Extraction of new solitary wave solutions in a generalized nonlinear Schrödinger equation comprising weak nonlocality

This article delves into examining exact soliton solutions within the context of the generalized nonlinear Schrödinger equation. It covers higher-order dispersion with higher order nonlinearity and a parameter associated with weak nonlocality. To tackle this equation, two reputable methods are harnessed: the sine-Gordon expansion method and the G′/(bG′+G+a)-expansion method. These methods are employed alongside suitable traveling wave transformation to yield novel, efficient single-wave soliton solutions for the governing model. To deepen our grasp of the equation’s physical significance, we utilize Wolfram Mathematica 12, a computational tool, to produce both 3D and 2D visual depictions. These graphical representations shed light on diverse facets of the equation’s dynamics, offering invaluable insights. Through the manipulation of parameter values, we achieve an array of solutions, encompassing kink-type, dark soliton, and solitary wave solutions. Our computational analysis affirms the effectiveness and versatility of our methods in tackling a wide spectrum of nonlinear challenges within the domains of mathematical science and engineering.


Introduction
Nonlinear phenomena manifesting in a fields such as engineering, physics, computational mathematics, chemistry, and biological sciences find precise representation and simulation through nonlinear evolution equations (NLEEs) [1].The importance of nonlinearity in wave dynamics is unmistakable within the realm of nonlinear sciences.Recent years have witnessed a concerted research endeavor aimed at unraveling analytical solutions, particularly solitary wave solutions, for NLEEs.To efficiently solve NLEEs and provide insights into intricate procedures, researchers have developed a variety of tools, including computational algorithms, analytical and numerical techniques.Several methods for analyzing these data have been proposed, such as modified generalized Riccati equation mapping approach [2], Lie symmetry analysis method [3], Hirota's method [4], inverse Fourier transform [5] and many others.Solitons, also known as invariant or solitary wave solutions, are stable, non-dispersive wave events that preserve their shape and speed as they propagate.These solutions, which are present in many physical systems, including nonlinear optics and water waves, are distinguished by their capacity to hold their original shape even when subjected to particular changes.Notably, solitons usually do not show singularities, which adds to their importance in comprehending and simulating wave behavior.Many methods have been developed to find soliton solutions.Researchers have used a variety of mathematical methods, such as Kumar et al. [6] employed generalized exponential rational function technique to obtain the closed-form soliton solutions, Rani et al. [7] used Lie symmetry analysis to attain explicit solutions, Dhiman and Kumar [8] utilized symmetry reduction to analyze the dynamics of invariant solition, exp (−ϕ(ξ)) expansion method is applied by Mathanaranjan [9] to find periodic and dark solitons, Zhao et al. [10] used exponential function technique to obtain soliton solutions and many others [11][12][13].
Erwin Schro ¨dinger introduced the Schro ¨dinger equation in 1925, which was formally published in 1926, laying the cornerstone for his subsequent scientific pursuits.The nonlinear Schro ¨dinger equation (NLSE) holds a prominent position among the remarkable nonlinear scenarios emerging from mathematical modeling.Numerous scientific domains, spanning mathematical finance, optical fibers [14], biology, fluid dynamics [15], quantum mechanics [16], plasma physics [17], quantum field theory [18], and the propagation of solitary waves in piezoelectric materials [19], to name a few, encompass significant applications for this equation.Diverse effective methodologies, including the Kudryashov method [20], the unified technique [21], the sine-Gordon expansion technique [22], the exponential method [23], and the extended trial equation method [24], have been extensively employed in the formulation of soliton solutions for nonlinear partial differential equations, with a particular emphasis on NLSEs.As cited in reference [25], an overarching NLSE encompasses both Kerr nonlinearity and fourth-order dispersion effects, exerting influence over the propagation of light waves.
Within this particular model, the variables t and x signify the temporal and spatial dimensions, respectively.Meanwhile, the symbol P is employed to denote the complex envelope of the electric field.The cubic nonlinearity coefficient is denoted by the parameter C 2 , and the coefficient accounting for fourth-order diffraction or dispersion is represented as C 1 .When dealing with the propagation of small wavelengths in optical transmission materials, it becomes notably important to factor in nonlocal effects.As elucidated in reference [25], the aforementioned model necessitates refinement through the incorporation of an additional component to account for the influence of weak nonlocality.
In this context, the term @ 4 P @x 4 characterizes the impact of fourth-order dispersion, while |P| 2 P encapsulates the Kerr nonlinearity.Furthermore, the term P(|P| 2 ) xx signifies the presence of weak nonlinearity.Recently, this model has been employed for the investigation of optical pulse localization within guided wave structures, particularly in scenarios involving a quasiperiodic linear factor [26].
Our motivation, therefore, centers on the quest for soliton solutions within the aforementioned models.We intend to achieve this goal by leveraging recently established and reliable techniques, such as the sine-Gordon expansion technique and the G 0 =ðbG 0 þ G þ aÞ-expansion method.These methods promise to address lingering questions from previous studies [27,28].The main aim of the sine-Gordon expansion method is to seek analytical solutions for a specific generalized nonlinear Schro ¨dinger equation.Its core objective is to simplify the equation while preserving its fundamental characteristics.This standardization process paves the way for the discovery of exact solutions, thereby furnishing valuable insights into the behavior of the system described by the PDE [29].In addition, we employ another method in this piece of research paper to obtain exact soliton solutions for a gNLSE through the G 0 =ðbG 0 þ G þ aÞ-expansion method.This method offers a wide spectrum of solutions, encompassing singular exponential, trigonometric, and kink soliton solutions, as detailed in reference [30].The discovered solitons have a crucial role in atmospheric research for controlling atmospheric gravity waves, maintaining wave stability, and preventing dispersion.Long-range coherent wave structure maintenance is essential for comprehending global atmospheric circulation patterns and how they affect weather occurrences.The effectiveness of the proposed methods in comparison to existing methods is determined by their problem-specific success and ability to provide accurate solutions, with factors such as ease of implementation, generality, and computational efficiency playing important roles in their evaluation.Consequently, our study employs mathematical computations and graphical analysis to scrutinize the physical characteristics associated with these waveform shapes [31].
The significance of this study, which includes weak nonlocality, lies in its ability to offer a robust mathematical framework that effectively captures the intricate interplay between nonlinearity and the subtle effects of weak nonlocality.This equation is particularly pertinent in the investigation of kink, dark, and novel solitary waves within nonlinear optical systems.Its value extends to serving as a valuable tool for comprehending and forecasting the behavior of solitary waves across diverse physical and optical scenarios.
The paper's structure unfolds as follows: In Section 2, we embark on a mathematical scrutiny of the model, aiming to derive the ordinary differential equation (ODE).Sections 3 and 4 delve into an extensive exploration and practical application of the sine-Gordon expansion method.Moving forward, Sections 5 and 6 shed light on the fundamental elements and practical utility of the G 0 =ðbG 0 þ G þ aÞ-expansion method.Section 7 presents a comprehensive overview of graphical illustrations, enriching our grasp of the subject matter.Lastly, Section 8 encapsulates our findings and contributions, serving as a concluding section that distills the pivotal insights of this study.Moreover, we also discuss future directions.

Section 2: Mathematical analysis
For crafting solutions within this specific model, we employ the following Ansatz: Swapping the expression from Eq (2) into Eq (1) generates the following ODE:

Section 3: Analysis of the sine-Gordon expansion method
Contemplating the sine-Gordon equation presented below: Here, we have the function denoted as P(x, t), and 'e' signifies any non-zero real number.
Step 1: Upon employing the traveling wave transformation P(x, t) = p(z), where z = x − Xt, Eq (4) simplifies into the subsequent nonlinear ODE.
While considering the function as P defined by p(z), where z represents the amplitude of the traveling wave, and X signifies the velocity of wave propagation, integration of the aforementioned equation results in the following expression: Here, B represents the constant of integration.Taking B = 0, p 2 ¼ WðzÞ and h 2 ¼ e 2 1À X 2 .Eq (6) assumes the subsequent form.
When we opt for h = 1, the equation above simplifies to the following form.
Applying the separation of variables method to solve Eq (8), we derive two significant relationships as follows: Step 2. The anticipated solution for Eq (5) can be expressed as follows.
pðzÞ ¼ Using Eqs ( 9) and (10), the above Eq (11) becomes To determine the positive integer ℵ, we utilize the balancing principle, a method that considers the balance between the highest-order derivative and the most influential nonlinear term found within the nonlinear ODE.
Step 4: Upon inserting Eq (12) into the transformed ODE, we arrive at an algebraic equation comprising various powers of sin(ϑ)cos(ϑ).By setting the coefficients of each power to zero, we establish a system of algebraic equations featuring the variables A 0 , A j , and B j .By solving this resulting system for A 0 , A j , and B j and subsequently substituting these values into Eq (11), we can construct a range of exact solutions for the transformed ODE.

Section 4: Applications of the sine-Gordon expansion method
In this section, we employ the sine-Gordon expansion method (SGEM) to solve Eq (3) and derive soliton solutions.Once we apply the balancing procedure, the solution takes the following form.
Employing the solution methodology as elucidated in Section 3, we calculate the values of the unknowns A 0 , A 1 , and B 1 , yielding the following results.
SET 1: By utilizing the values of A 0 , A 1 , B 1 , β, and α obtained from SET 1 in Eq (13), and subsequently substituting them into Eq (2), we observe the solution as follows: SET 2: Similarly, when we employ the values of A 0 , A 1 , B 1 , β, and α from SET 2 in Eq (13), and subsequently substitute them into Eq (2), the solution appears as follows.
Section 5:The G 0 =ðbG 0 þG þaÞ-expansion method Suppose that the transformed ODE possesses a solution in the following manner.
Here, we have V ¼ VðzÞ ¼ G 0 bG 0 þGþa , with a and b being nonzero parameters.The value of ℵ can be ascertained using the principle of homogeneous balancing.Additionally, we introduce arbitrary constants denoted as F j , which will be determined subsequently.Furthermore, V = V(z) satisfies the subsequent ODE.
Additionally, we have G = G(z) which represents the solution of the following ODE.
In Eq (15), we encounter real numbers λ and μ.This equation exhibits the following categories of solutions, Type 1: ðÀ lÀ ffi ffi ffi

:
We can alternatively express V = V(z) as follows: By inserting Eqs ( 14) and ( 15) into the transformed ODE and setting the coefficients of V j to zero, we establish a system of equations.Solving this resultant system of equations enables us to ascertain the values of the arbitrary constants.

aÞ-expansion method
In this subsection, we employ the G 0 =ðbG 0 þ G þ aÞ-expansion method to extract soliton solutions from Eq (14).Following this method, for the case when N = 1, the assumed solution for Eq (14) takes the subsequent form.
Here, we have the unknown constants represented as F 0 and F 1 .Employing the G 0 =ðbG 0 þ G þ aÞ-expansion method results in the following set of solutions.SET 3 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi The solutions corresponding to SET 3 are presented below.In accordance with Type 1, when Θ = λ 2 − 4μ > 0, we obtain the following results 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where, (18) becomes 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where, 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where, In alignment with Type 2, when Θ = λ 2 − 4μ < 0, we observe the following outcomes.
2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where ffi ffi ffi ffi ffi ffi ffi ffi À Y p r 1 ¼ 0, then Eq (21) becomes 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where Moreover, if ðl À 2Þr 1 À ffi ffi ffi ffi ffi ffi ffi ffi À Y p r 2 ¼ 0, then Eq (21) becomes 2 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where Section 7: Graphical illustration Graphs serve as powerful tools for conveying information and effectively illustrating the solutions to the problem in a logical manner.In this context, a series of graphs representing solutions for various parameter values were presented by using the SGEM and the G 0 =ðbG 0 þ G þ aÞ-expansion method.The SGEM initially acts as a valuable instrument for revealing exact solutions across a wider range of nonlinear Schro ¨dinger equations.Subsequently, we employ Mathematica 12 software to create visual representations of these exact solutions, specifically for both SET 1 and 2. Within these solutions, established through the SGEM, we encounter hyperbolic functions that encompass real and complex values.The outcomes identified in our analysis can be harnessed to assess real-world characteristics related to wave amplitudes and widths via these hyperbolic function solutions.The hyperbolic function solutions offer valuable insights into critical parameters such as wavelength and frequency, as exemplified in the equations provided earlier.It's worth noting that as wavelengths extend, their impact on a global scale becomes increasingly significant.Consequently, delving into the mathematical foundations of these natural phenomena becomes imperative, whether to mitigate their potentially destructive effects or harness them as valuable energy sources.
Through the utilization of 3-D and 2-D graphs, this study effectively demonstrates the behavior of the real and imaginary components of the computed solutions, designated as P 1 (x, t) and P 2 (x, t).Specifically, the figures labeled as Figs 1 and 2 present the contextual values of the specified parameters: A 0 = B 1 = 0, C 2 = 1, C 3 = 4, and C 1 = −1.Our study yields a diverse array of solutions, including Dark, Kink, anti-Kink, and single Soliton solutions.In the realm of optics, NLSE plays a pivotal role by elucidating the intricate behavior of optical pulses within nonlinear media.Soliton solutions, with their unique characteristics, particularly for nonlocal Schro ¨dinger equations in optics, underscore the significance of our findings.
The application of the G 0 =ðbG 0 þ G þ aÞ-expansion method in our study proves to be an effective approach for discovering solitary wave solutions within a more extensive range of NLSEs.The graphical depictions notably illustrate these solutions' solitary waveforms.Upon employing specific parameters, such as λ = 2.2, μ = 1, r 1 = r 2 = 1, C 1 = C 2 = 1, and C 3 = −1, within our solution P 3 (x, t), as exemplified in Fig 3, novel soliton solutions are generated.Likewise, when we suggest values like λ = 1.2, μ = 1, r 1 = r 2 = 1, C 1 = 1, C 2 = 1, and C 3 = 1 in our solution P 6 (x, t), illustrated in Fig 4, we obtain new solitary wave solutions.Utilizing Mathematica 12, we created graphical representations of the exact solutions mentioned above.These graphs encompass both real and complex solutions derived from the new soliton solution and are presented in both 3D and 2D formats.

Section 8: Conclusion
This article focuses on exploring exact soliton solutions for the generalized nonlinear Schro ¨dinger equation.The equation was characterized by higher-order dispersion with higher order nonlinearity, and a parameter associated with weak nonlocality.To address this challenging equation, we employed two well-established techniques: the sine-Gordon expansion method and the G 0 =ðbG 0 þ G þ aÞ-expansion method.By applying these methods in conjunction with appropriate traveling wave transformations, we successfully derived novel and highly efficient solitary-wave solutions for this model.To gain insights into the equation's physical implications, we harnessed Wolfram Mathematica 12 for generating both 3D and 2D graphical representations.These visualizations offer profound insights into various aspects of the equation's behavior.Through manipulation of parameter values, we unveiled a variety of solutions, including kink-type, dark soliton, and solitary wave solutions.Our computational analyses underscore the effectiveness and versatility of these methods in addressing a wide range of nonlinear problems in the domains of mathematical science and engineering.Sometimes these techniques, may not be as reliable in producing accurate solutions due to their limited applicability to a diversity of equations and possible convergence problems.Future prospects involve the exploration of alternative analytical techniques, including the derivation of   generalized solutions like Jacobi elliptic function solutions and multiple solutions.Furthermore, conducting bifurcation analysis, chaos analysis, and investigating modulation instability can be pursued.This research may lead to the discovery of N-soliton solutions, rogue waves, breathers, and various other intriguing solutions.